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Abstract 

Normal behavior of the thermodynamic properties of a Fermi gas in d > 2 dimen- 
sions, integer or not, means monotonically increasing or decreasing of its specific heat, 
chemical potential or isothermal sound velocity, all as functions of temperature. How- 
ever, for < d < 2 dimensions these properties develop a "hump" (or "trough") which 
increases (or deepens) as d — > 0. Though not the phase transition signaled by the 
sharp features ("cusp" or "jump") in those properties for the ideal Bose gas in d > 2 
(known as the Bose-Einstein condensation), it is nevertheless an intriguing structural 
anomaly which we exhibit in detail. 
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1 Introduction 

The familiar ideal Fermi gas is revisited for any positive space dimension, d > 0. Ideal 
Fermi gases (IFG) have been discussed thoroughly in general by several authors [l]-[4] and a 
detailed study of the quantum behavior in any dimension at sufficiently low temperatures in 
these systems have also gained interest as possible precursors of a paired-fermion condensate 
at lower temperatures, these have been studied experimentally in ultra-cold fermionic clouds, 
e.g., with fgK neutral atoms in opto-magnetic traps [5]-[8]. From the thermodynamic grand 
potential we generalize the system to any dimension, we calculated the thermodynamics 
properties and analyzed the results for the chemical potential, the heat capacity and the 
isothermal velocity of sound where we found that there are structure if li < 2 in contrast 
with the ideal Boson gas where there is structure if o? > 2 [9, 10]. In Sec. 2, we calculate the 
thermodynamic grand potential for the non-interacting Fermi gas in d-dimensions, we find 
the thermodynamic properties of these systems and deduce a generalized density of states 
(DOS). In Sec. 3, we calculate the chemical potential /i(T), the internal energy and the 
specific heat as function of absolute temperature T for any positive space dimension d and 
analyzed the results obtained for c? < 2. In Sec. 4 we obtain the isothermal and adiabatic 
velocities of sound for these systems. Sec. 5 contains our conclusions. 



We consider an ideal quantum gas of fermions in d > dimensions of mass m in vacuum with 



2 Ideal Fermi gas in d > dimensions 



a quadratic dispersion relation. The Hamiltonian is then H — X]^=iPf/2m, its eigenvalues 
are given by 




(1) 



j=i 
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where L is the size of the "box", n = (ni, n2, nd), and where Ui — 0, ±1, ±2, .... Since 
ki = {27r/ L)ni, (1) can then be rewritten as 

fc2 d 
1=1 

The thermodynamical properties of this system follow from the thermodynamic (or grand 
potential) ^1(7", L'^, ji) — U — TS + jiN, where U is the internal energy, T the absolute 
temperature, S the entropy, ii the chemical potential, and N the number of particles. We 
may write 0(T, L*^, fi) in generalized form as (see p. 134 of [1]) 

Q(T,L^//) = -A;bT J]ln[l + e-^(^'=i-'^)], (3) 

with (5 = l/ksT, and is the Boltzmann constant. Using the logarithm expansion ln(l + 
x) — — S^i(— 2^)V^ valid for x < 1, (3) becomes 



ki 1=1 

CO 



I 



1=1 ki 



In the continuous limit where ti?/mLP' <^ ksT, the summation over ki can be approxi- 
mated by integral, namely J2ki — (2s + l)(L/27r)'^/ d'^/cj. Thus 



n{T,V,fj,) = A;Br(2s + l)(L/27r)'^E ^ ^ / e"'^'^'^"/^'")^" x 

00 J —00 

with s the particle spin. Integrating, Eq. (5) becomes 

The infinite sum can be expressed in terms of the Fermi functions fa{z) (see Appendix D of 
Ref. [1]), 



U{z) = -Y = ^ / dx—^ , (7) 



so that 

d/2 



which defines A^, and where z = e^^^^\ 



Prom (8) it is possible to find the thermodynamic properties of a monatomic gas using 
the relation dO, — —SdT — PdV — Ndji, where V = L*^ is the system volume and P its 
pressure. In this representation, the grand potential $1(7, y,/i) is the fundamental relation 
leading to all the thermodynamic variables of the system, namely, 

3 Thermodynamic variables 
3.1 Chemical potential 

We determine the chemical potential /u(T) as a function of absolute temperature T for any 
positive space dimension d > from the number equation which is given from (8) and (9), 



^ = -^(|/../«w)^,, = ^.^. (10) 

which after substituting Aa in (10) and using (7) becomes 



r{d/2) \2^) X 60(^-1^) + 1 ^^^^ 

Since N — deg{e)n{e) where g{e) is the density of states (DOS) and n{e) = [e^^^"'') + l] ^is 
the Fermi-Dirac distribution, then 

Since [e^{e-M(T)} + ^j-^ Qi^Ep-e), with ^(x) the unit step function, /i (0) = Ep = 

t^k'^p/2m the Fermi energy, kp being the Fermi wavenumber, we see from (11) that 
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and we recover the expressions obtained in [3, 11] for the fermion number density with 
s = 1/2, i.e., 

_ N _ kf 

2'^-27r^/2rfr(rf/2)' ^ ^ 

which gives the famihar results 2/^^/71, kp/27i and kp/3n'^ for d = 1, 2 and 3, respectively. 
Defining Ep = ksTp, the chemical potential /x(T') is obtained from (11) and (13) by solving 
numerically the following equation 

T^I^V{d/2)U,2{e^^) = {2/d) T//\ (15) 

For d — 2, fi{z) from (7) is the ordinary log function so that (15) gives 

/i(^)=ln(l + ^)=TF/T, (16) 

leading to the relatively well-known explicit [12] formula 

l,(T)/Ep^T/Tp\n{e^-/^-l) ^ 1, (17) 

which is clearly not expandable in powers of T/Tp as in the so-called Sommerfeld expansion 
[13]. However, for d ^ 2, ijl{T) is not an exphcit, closed expression in T/Tp; numerical 
analysis is required to extract it. We now determine the chemical potential ^x{T) as a 
function of absolute temperature T for any positive space dimension, d > 0. This requires 
solving (15) numerically when d ^ 2. Whereas /i(T) in both = 2, 3 is well-known to 
decrease in T from the constant Fermi energy Ep (which depends only on the number density 
of fermions and on their mass) at T = 0, towards the well-known classical value diverging 
logarithmically to — oo, for all d < 2 we find novel, anomalous behavior consisting in a 
curious temperature non-monotonicity: /i(T) first increases quadratically as T is increased, 
then changes curvature, acquires a maximum, and finally decreases monotonically to the 
classical value. This translates into a "shoulder" or "hump" in the heat capacity as function 
of temperature, and may be relevant in the study of quantum "dots," "wires" and "wells" 
of modern-day opto-electronics [14, 15]. 

Results are exhibited in Fig. 1 for several values of d. The unexpected rise in fi{T/Tp) 
for c? < 2 with increasing T/Tp is novel and anomalous, though it was reported in Ref. [3] 



and more completely in Ref. [4] albeit with some errors. For c? = 1 it is graphed in Ref. 
[16] p. 192 for low temperatures, but without comment. Using the large- 2; expansion (Ref. 
[1] p. 510) f^(z) ~ (\nzy/r(a + 1) + (n'^/G) (a - 1) (ln^)^~^ /r(a) + • • • the d-dimensional 
Sommerfeld expansion for d 7^ 2 becomes 

t,{T)/E^ --^ l-{d/2-l){ny6){T/T^f + 0{T% (18) 

Note that the first correction to unity is positive for all ci < 2, and since for large enough 
T the chemical potential must diverge negatively to approach the classical value a "hump" 
will emerge. Using the small-z expansion fa{z) ~ 2; — H , one gets 

H{T)/E^ -{T/T^)\n\v{d/s){T/T^f\d/2)\, (19) 

or the well-known classical limit for large T. 

Prom (8) and (9) we obtain this expression for the entropy S in dimensionless form 

S/Nks = {d/2 + 1)1^1^ . d^(^) . 0, (20) 

where the T ^ limit is obtained by using the large- 2; expansion for fd/2{z), and clearly 
complies with the third law of thermodynamics. 

3.2 Internal energy 

The internal energy can be obtained from (see p. 159 of Ref. [1]) 

t,(r.v,^-*.r^A(_^)^^^. pi, 

Substituting (8) in (21) and comparing with (8) we find that (recaUing that V = L'^) 

d Aa 

2/?^ 



U(T,V) = '-^J,/,,,(z)^'-n^'-PV. (22) 



Eq. (22) is a generalization of the relation PV = |t/ for an ideal gas of fermions (and in 
fact, also bosons), in the nonrelativistic limit. Comparing (10) with the last equation we 
find 



For T ^ 0, using the large-z expansion for fd/2{z) one has 

In the classical limit T — >• oo, 2; << 1, we again use fa{z) ~ 2; — + ■ ■ ■ and from (23) 

get 

U/N—^ %T, (25) 

1 — >oo Z 

in accordance with the equipartition theorem in d dimensions. 



3.3 Specific heat 

The specific heat at constant volume Cy = [dU{T, V)/dT]j^y follows from (23) 



Substituting (10) in the last equation, this becomes 
where we used the relation 



(26) 



(27) 



-[w) --^Bm^)P^y (28) 
^ \^ ) Ny fd/2-i{z) 

which is extracted from the (vanishing) derivative with respect to T of the number equation 
(10). 

Equations (24) and (25) for large-and small-^ may be differentiated to yield the specific 
heat as 

Cy/Nks ^ (^) (29) 



and 



2Cy/dNkB^ l + ^{l-d/2). (30) 



Fig. 2 (top panel) illustrates the "hump" developed by the Cy (T) for all c? < 2. These 
peculiar results for the ideal Fermi gas in 0? < 2 dimensions, manifesting "structure" in the 
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form of anomalous (non-monotonic) behavior in the chemical potential and in the specific 
heat, contrasts sharply with the ideal Bose gas where the "structure" appears for all d > 2, 
(see Fig. 2.5 of Ref. [9] for integer d, and Fig. 2 of Ref. [10] for all d > 2). The (sharp) 
structure observed is the Bose-Einstein condensation (BEC) whose signature is a "cusp" in 
the specific heat at the critical transition temperature for all 2 < d < 4, and a "jump" in its 
value there for all 4 < o? < oo. 



4 Sound velocities 

To further exhibit the anomalous behavior for ci < 2 we have also calculated mechanical 
properties, e.g., the isothermal ct and adiabatic cs velocity of sound in the IFG. These are 
defined as 

where from (8) the pressure P = —D,/V is given by 

P = iA./L-') ^^pW. (32) 

4.1 Isothermal sound velocity 

From (31) and (32) one has 

'"4m = ^(^)^. 

Using the relation 

which is obtained from (10), we find 

mcUT) = ksTJ^^. (35) 

Jd/2-l[Z) 
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Normalizing with Ep = mvp/'i = ksTp (35) we have 



,2 1/^/^ X fd/2{z) 



[ct{T)/vpY = -{T/Tp)-P^ (36) 
Using the large- 2; expansion for fd/2{z) as T — > we find 

[cT{T)/vpr /^/ (c^^^) + ^ - ^^1^ (37) 



^ ^ 12 V 2 J n/E^ 
or, since at T/Tp — 0, /i — Ep one gets the famihar result 

ct(0) = vp/y/d. (38) 

In the limit of small fa{z) — > 2; as T — > 00 we have 

[ct{T)/vp]\—^ T/2Tp, (39) 

which corresponds to the classical limit. Prom (37) and (39), we conclude that ct{T)/vf 
decrease in T for all d < 2 from 1 at T = 0, then changes curvature, acquires a minimum, 
and finally increases hnearly as T to the classical value. In Fig. 3 we plot Cy(T)/cf.(0) 
where Ct{T) is the isothermal sound velocity for an IFG as function of temperature T jTp 
for c? = 1/4, 1/2, 1, 2 and 3. Note that ct develops a "trough" for c? < 2 which deepens as d 
decreases, more clearly seen in the bottom left panel where we show the minimum values in 
the isothermal sound velocity compared with ct(0) as function of Log d, and on the values 
of the temperature Tmin corresponding to the minimum in top panel. Thus c|.(T) never 
becomes negative so that ct(T) is always real. 

4.2 Adiabatic sound velocity 

The adiabatic sound velocity cg in a Fermi gas is obtained from the adiabatic state equation 
for an ideal gas in d-dimensions, namely (Ref. [1], p. 229) 

PV^ = const = TV^c, (40) 
where 7 = l-|-2/d, A'"is the particle number and c some constant. Hence 

P = v?c. (41) 
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Using (41) in the second equation of (31) we find 



^1 - (-] 
^ \dn) 



N,S 



n 



(42) 



Substituting 7 and (41) in (42) we fiave 




( 



d J n 




( 




fd/2+l{z) 
fd/2{z) 



(43) 



from which is obviously monotonic for all d. 

5 Conclusions 

After constructing the grand potential, thermodynamic properties were determined along 
with densities-of-states for an ideal Fermi gas in d dimensions, integer or not. For < c? < 
2 dimensions these properties develop a hump in the chemical potential /Li(T) and in the 
specific heats Cv{T) and a trough in the isothermal velocities of sound ct{T). This structure 
contrasts with the case of the ideal Bose gas structure ocurring, however, for d > 2 commonly 
associated with BEC, and characterized by a cusp in the specific heat when 2 < d < 4 and 
a jump for all 4 < 0? < 00. 
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Figure 1: Chemical potential (in units of Ep) for an IFG in o? = 1/2, 3/4, 1, 2 and 3 

dimensions as functions of temperature T (in units of Tp). On the left bottom panel we 

show the maximum values in the chemical potential as function of d, and on the right the 

values of the temperature T^ax corresponding to the maxima in top panel. 
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Figure 2: Specific heat for an IFG in d — 1/4, 1/2, 1, 2 and 3 dimensions. On the left 
bottom panel we show the maximum values of the specific heat as function of d and on the 
right the values of the temperature Tmax corresponding to the maxima in top panel. 
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Figure 3: The quantity c^{T)/c^{Q) where ct{T) is the isothermal sound velocity for an IFG 

as function of temperature T/Tp. On the left bottom panel we show the minimum values in 

the isothermal sound velocity compared with ct(0) as function of Log d, and on the values 

of the temperature T^in corresponding to the minimum in top panel. 
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